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Abstract: In this paper, we study central and non-central limit theorems for partial sum of 
functionals of general stationary Gaussian fields. We apply our result to study drift parameter 
estimation problems for some stochastic differential equations related to stationary Gaussian pro¬ 
cesses. 
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1 Introduction 

While the statistical inference of Ito type diffusions has a long history, the statistical analysis for 
equations driven by fractional Brownian motion (fBm) is obviously more recent. The development 
of stochastic calculus with respect to the fBm allowed to study such models. We will recall several 
approaches to estimate the parameters in fractional models but we mention that the below list is 
not exhaustive: 

• The MLE approach in on, m- In general the techniques used to construct maximum 

likelihood estimators (MLE) for the drift parameter are based on Girsanov transforms for 
fBm and depend on the properties of the deterministic fractional operators (determined by 
the Hurst parameter) related to the fBm. In general, the MLE is not easily computable. 

• A least squares approach has been proposed in M- The study of the asymptotic properties of 

the estimator is based on certain criteria formulated in terms of the Malliavin calculus (see 
m- In the ergodic case, the statistical inference for several fractional Ornstein-Uhlenbeck 
(fOU) models has been recently developed in the papers [Lf], El, [3] , ff3| . [To] , [9]. The 
case of non-ergodic fOU process of the first kind and of the second kind can be found in (4j 
and m respectively. 

Our aim is to bring a new idea to develop the statistical inference for stochastic differential 
equations related to stationary Gaussian processes by proposing a suitable criteria. Our approach 
is based on Malliavin calculus and it makes in principle our estimator easier to be simulated. 
Moreover, the models studied in na, ei, 0 , mi become particular cases in our approach (see 
section 4). 
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2 Elements of Gaussian analysis and Malliavin calculus 

For the essential elements of Gaussian analysis and Malliavin calculus that are used in this paper 
see the references m, Eg). 

Now recall that, if X, Y are two real-valued random variables, then the total variation distance 
between the law of X and the law of Y is given by 

d TV (X, Y) = sup | P [X € A] - P [Y € A] |. 

AeB(R) 

If X, Y are two real-valued integrable random variables, then the Wasserstein distance between 
the law of X and the law of Y is given by 

d W {X, Y) = sup \Ef(X)-Ef(Y)\ 

f&Lip{ 1) 

where Lip( 1) indicates the collection of all Lipschitz functions with Lipschitz constant ^ 1. 

The following well-known direct consequence of the Borel-Cantelli Lemma (see e.g. [IS]), will 
allows us to turn convergence rates in the p-th mean into pathwise convergence rates. 

Lemma 1 Let 7 > 0 and po € N. Moreover let (Z n ) n& ^ be a sequence of random variables. If 
for every p ^ po there exists a constant c p > 0 such that for all n E N, 

(E| Z n \P) 1/p ^ c p ■ n~\ 

then for all e > 0 there exists a random variable i] £ such that 

| Z n \ ^ r] £ ■ n~ 7+£ almost surely 
for all n € N. Moreover, E|%| p < 00 for all p ^ 1. 


3 General context 


3.1 Central and non-central limit central theorems 

Consider a centered stationary Gaussian process Z = (Zfc) fcgZ with covariance 

rz{k) := E(ZoZk) and \rz(k)\ ^ 1 for every k € Z. 

Let us define the following normalized process 


We also define 


Y k : = 


Zk 

y/rz{ 0)’ 


k € Z. 


q/2 

V q , n (Z) = k,n) 

k=1 


(1) 

( 2 ) 


2 




where 


n— 1 

f2k,n ■= d g , 2k (Z)^J2 £ T 2k 

with Yi = Y(si) and d q , 2 k(Z) are constants. 

If |ry(fc)| 2 < oo, we dehne 

v q (Z) : = lim E [V q 2 n (Z)\ . 

n —^oo L J 

We have the following central and non-central limit theorems. 


Theorem 2 Let (V 2 , n (Z)) n>0 be the sequence defined in {2j), we let Y^ = 
covariance r Y {k) = E(YqY}L). Denote N A7(0,l). 

1) Then there exists a constant C > 0 depending on q and rz( 0) such that, 


Zk 

y/rz(0) ’ 


k € Z with 


dry 


V q , n (Z) 


E [V q 2 n (Z)\ 


■-,n\ ^ 


c 


E 


VU Z ) 


E 


Vi AY) \ ^(V 2 ,n(Y)) + KfiV 2 ,n(Y)). (3) 


We also have in the case when, 'Zke'Z l r r(^’)| 2 < °°> 


dry 


Vg,n(Z) 

VMZ)' 


N 




C 


v q (Z) 


r~ 


r 

vi,AY) 


Vln(Y) jK 4 (V 2 , n (Y)) + K4(V 2 , n (Y)) + \V g (Z) - E [V^Y )] | 


(4) 


In addition: 

• VY.kezAri.k)) 2 < oo, 


k*(V 2 , n (Y)) < Cn~ l ( Y, (r Y (k)f 3 

Jfc|<ra 


Ifry(k) = k 2 , 


K4{V 2 ,n(Y)) 




C 


V 2 2 n (Y) 


log (n)’ 


where E [V£ n (Z)] ~ 4 d? q 2 {Z) log(n). 

2) Assume that r Y (k ) = \k 

V q , n {Z) law^ d q , 2 (Z) 

n h~ a Ad 


with 0 < a < ^. Then 


F n 


(5) 


( 6 ) 


(7) 


3 




























where D = e iy \y\ a 1 dy = 2r(a) cos(^), and 


F — 
1 oo — 


A(x+y) _ i 

i{x+ y ) - \xy\^W(dx)W(dy) 


i{x + y) 

with W is a random spectral measure of the white noise on R. 


Proof. Since 


Vq,n{Z) 

^E[V£ n {Z)\ 


1 ) 1.2 


, t»y PH Proposition 2.4] we have 


drv 


V q ,n( z ) 

V E R»( Z >]’ 


N \ < 2E 


1 -ID 


V q , n (Z) 


-DL 


-i V q , n (Z) 


yjE [V q 2 n (Z)] ’ y/E [V£ n (Z)] 


H 


On the other hand, exploiting the fact that 


E 


E 


(• hk{f2k,n )) — (Dl2k(f2k,n),—DL I 2 k{f2k,n))H 


= 0 


we obtain 

E 


1 - >n fkB . DL- 




[V£ n ( Z)] ’ ^ KniZ)] 

1 / ?/2 /- 

r -—-T ^ ((2fc)- 1 || J DI 2fe (/ 2M )||^) + Y, (20“ 1 ^K^2fc(/2fc,n). £>/ 2l(/2t > n))w| 

V fc=1 l<fc<Z< 0 /2 


+ £ (20- 1 ^|(^/2fc(/ 2 M.).^(/2/,n))«| • 

l^l<k^q/2 ) 

Moreover, by pa Lemma 3.1] we have 

_ 2fc—l 

yJVar((2k)-i\\DI 2k (f2k,n)\\ 2 n ) = (2k)~ 2 Y ft? 

3 = 1 

and for k < l 


(4 A) 2j)!||/ 2 fc jTl (i)/ 2 fc jn || wg i4 fc _2 3 -, 
i 


E 


{(2l)-\DI 2k {f 2kin ),DI 2l (f 2lin )) H ) 2 } < (2A;)! fc) {2l-2k)\E\{I 2k {f 2Kn )f]\\f 2Kn ® f 2k 


21—2k 


2k—1 


+2k 2 — l)- 2 (j-l 1 ) (i-! 1 ) "b 2/ — 2j)! ^ ||/ 2 fc,n 2 ® ./2fc,n||%®2j + ||/2fc,n^p ./2fc,rillft®2j 


J=1 


2l-j 


Combining this together with the fact that for every 1 ^ s ^ 2k — 1 with k € {1,..., q/2} 

n 

||/2fe,n®/2fe,n||^®4fc-2 S ^ d q2k {Z)n~ 2 Y r Y( k i - k 2 )r^(k 3 - k 4 )ry ~ s (ki - k 3 )ry~ s (k 2 - fc 4 ) 

^1,^2,^3,^4 = 1 

n 

< dq t2k (Z)n ~ 2 ^ r Y (h - k 2 )r Y {k 3 - k 4 )r Y (ki - k 3 )r Y {k 2 - k±) 

kl ,^2 ,&3 5^4 = 1 

= <2fc(^)«4(^,n(^)) 


2 

%®4fc 
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we deduce that there exists a constant C depending on q and rz( 0) such that 
y/Var {(2k)-mDI 2k (f 2kjn )\\^) ^ Cy/K A (V 2 , n (Y)) 


and 


{2l)- 1 E\{DI 2k {f 2Kn ),Dh l {f 2ltn )) H \ < (E \{{2l)-\DI 2k (f 2Kn ),DI 2l {f 2l , n )) H ) 


1/2 


< C\ E 


(hk(f 2 k,n)r V K4(y 2l n(Y)) + K 4 (V 2 , n (Y)) 


^ c 




n—1 


n 1 Y2 r Y (* “ j)\J «4(V2,n(- Y)) + K4(V 2 , n (T)) 


i,j =0 


< C 


N 


n—1 


n 1 ^2 r Y (* “ j)\J Ki{V 2 ,n{Y)) + 


*,3=0 


Furthermore, 


dry 




E [V q 2 n (Z)\ 


■-,N\ ^ 


c 


n—1 


E [V q 2 n (Z)} \ 

c 

E [V 2 n {Z)\ 


n 1 y~] - j)^/K 4 (V2,n(^)) + K 4 (V2,n(^)) 


*,3=0 




<7 


E 


V?,n( Z ) 


E 


E 


V 2 f t n(X) J^(V 2 , n (Y)) + K A (V 2 , n (Y)) 


V? >n (Y) J^n(Y)) + «4 (V 2 , n (Y)). 


Thus the estimate © is obtained. By a similar argument we obtain ©. 
Moreover, if ( r y(^)) 2 < oo, it follows from [6] that 


«4 {V2,n{Y)) < 


Cn 


-l 


E 


vU Y ) 


£ (’vM) 4 '' 3 

|fc|<n 

3 




Cn 1 ( ^ (ry(fe)) 4/3 ] 


. |fe|<n 
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and if r Y (k) = k 2 for k E Z*, it follows from (l‘23[) that 

n— 1 9/2 


n—1 


E i V q,n( Z )\ = d qA Z )“ 2 M* “ j)| 2 + £ d q, 2 k( Z )^- E ^ Y (* “ ^)l 

i,j=0 k=2 i,j =0 


2k 


2k 


9 n_1 \ ra-i 

2 d^ 2 ( Z ) j 1 + ~E( n ~-?)-7 1 ) + '^2 d q,2k( Z )—^— E \ r Y(i~j) | 

j=l / k=2 i,j =0 

n-i n-1 \ 9/ 2 / 0 i\, «-l 

2d q,2( Z ) ( 1 + 2 E J 2 ) + 5Z d 9,2fe( Z )—~ X] l r H*-j)l 

j=l j=l / fc=2 ij'=0 


2fc 


~ 4d 2 2(Z)log(n). 
Thus, as in [7] we conclude that 
K4(V2,nP0) _ K 4 (V 2) nP0) 


( 8 ) 




^ A (0)^ ^n(^) 

1 


2 I T * 






r 


C 


2 E ry(fci - k 2 )r Y (k 3 - k 4 )r Y (k 1 - k 3 )r Y {k 2 - fc 4 ) 

>^2 >^3 5^4=1 


n 2 log"(n) 
C 


\0)E [F 2 2 n (Z) 

n 

E r H fc i “ k 2 )ry(h - ki)r Y (ki - k 3 )r Y {k 2 - fc 4 ) 


ki ,k2 ,k$ ,k4=l 


log (n) 


(9) 


Hence the inequalities © and © are satisfied. 

The convergence © is a direct consequence of [ III Theorem 1], 


3.2 Hermite variation 

Let g € N* be even, and let H q be the pth Hermite polynomial. So, H q has the following 


decomposition 


Define for every q ^ 2 even, 


#,(*) = E 




fc =0 


k\{q-2k)\2 k ' 


n —1 


Qq,n(Z) '■= ~ E Hq( Z k) 

n z —' 


( 10 ) 


fc =0 


and 

7z(g) := E [H q (Zo)] ■ 
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Then, we can write 


/ N_\ F (ryq-2k 

7Z 9 ’ E k\(q-2k)\2 kE \ Z ° 


ql 


E 


(-1)* 


2i / 2 ^ ife!(§ - fc)! 

fc=0 V2 ' 




2M§-* 


9 ! 


(t;(z 0 2 )-i) 9/2 , 


(f)!29/ 2 


We have the following almost sure convergence. 

Theorem 3 Suppose that Z is ergodic. Then, as n ^ 00 

Qq,n(Z) —■> 7 z(q) 

almost surely. 

Let Y % be the process defined in (JTJ). We have, 

9/2 

H q (Zi) - EHg(Zi) = J2 b q ,2 k H 2k (Yi) 


k=i 


where for any k E {1,... , % — 1} 


b q , q - 2 k(Z) = (-l) fc (rf(0)-r| ^0)) a q q _ 2 a q _l... a q q _f k +2 

l q-2k 


+ (-l) fc 1 (^1(0) - r| 2 (0)) a q q _ A a q q _l... a q q _f +2 

+ ... 

+ (-1) 1 (E(°) 


l q—2k 


and bq^g(Z) = 7^(0), with for every p even 

p!(—l) fc 


-,P — 


P~ 2k k\(p - 2k)\2 k 


k = 0,... ,p/2 


which verify 


p/2 


h p( x ) = E°p-2fc xP 2k - 


k =0 


Dehne 


( 11 ) 


( 12 ) 


V Qq jZ) := V^(Q,,n(^)-7z(9))- 
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As consequence, we can write 


where YJ = Y(si). 
Furthermore, 


Vq 9 JZ) = V^{Q q ,n{Z)-lz{q)) 


n— 1 


— YslHgiZi,) - EH q (Zi)] 

vn t=0 
-j n— 1 g/2 

b q ,2k{ Z ) H 2k(Yi ) 

v j=0 fc=l 
g/2 n—1 

fc=l V i=0 

?/2 / j n— 1 

J^ J 2fc fyf,2 k( Z )-r=J2 
k =1 V vn i=0 


r / <S>2fc 
22fc I 


E 



9/2 (ob\\ n ~^ 

E^EM-# 

fc=l *,j=o 

9/2 / 2 n—1 

E *?,»( z )( 2 *) ! 1 + - E<" - 1 - j)I^'0')I 


2k 


k= 1 
9/2 


i=i 

n—1 


n—1 


E'E»(Z)(2«! l + 2El>'y0)| a --E.'l’-l'0)l 


2 k 


k =1 


9=1 


i=i 


(13) 


We will need the following technical lemma. 


Lemma 4 Let (Zk)k^o be a stationary Gaussian sequence with E(Zq) < 00 , and let A > 0 and 
q £ N* even. Consider the sequence 


RqA A, Z fc ) := Lf g (Z fc - e- Afc Z 0 ) - LT g (Z fc ). 


Then for every p ^ 1 t/iere exits a constant c(X,q ) depending on \,q and E(Zq) such that 


Moreover for every e > 0 


j n—1 

~^2 R Q,q( X ’ Z k) 

1 k =0 


Lp(U) 


c(A,g) 

n 


1 

n £ 


n—1 


^RQ^Zk) 

k =0 


0 


almost surely as n 00 . 


(14) 


(15) 
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Proof. We have 


q—2i 


R Q .,(\,z k ) = E m ’ !( _ g, 2 , 

Combining this with the fact that Z is stationary and Gaussian, we obtain 


\\RQ,q(X, ^fe)lliP(f2) ^ c 0('V ( l) e 

Thus (fTlj) is obtained. 

For the convergence (I15[) , it is a direct consequence of (j!4[) and Lemma |T] ■ 
Applying Theorem [2] and Lemma 0] we conclude the following result. 


Theorem 5 Let (Lq 9 n (-^)) n>0 and (Rq^(X, Z k )) n>0 be the sequences defined in the above. 
1) Then there exists C depending on X, q and rz{ 0) such that 


dw 


r 

n 

I e 

V U Z \ 


( Q q ,n (Z k - e~ Xk Zo) - 7z (g)) , N 


< 


C 


E 


V lJ Z ) 


/ 

E 

%J Y \ 

\ 

A 

1 J TP 


^4(V Q2 JY)) + k 4 (V Q2 JY)) 

\ 


n 

+ v 


(.16) 


On the other hand if |ry (/c)| 2 < oo, we can write 

(< Q q ,n (z k -e- Xk z 0 ) -7z(g)) ,n) 


dw 


v Qq (Z) 


€ 


C 


v Qq (Z) 


V -^Z1 + X IE 

n 


v $ 2 JY ) sJ^(V Q2 JY)) + n 4 (V Q2 jY))+ v Qq (Z) — E 


Vi 


Qq,r 


where vq q (Z) := lim n _ ) . 00 .E Vq (Z) . Moreover, k 4 (VQ 2n (Y)) verifies (ED and /f6]j . 
2) Assume that ry(k) = \k\~ a with 0 < a < Then 


n 


( Qq,n (z k - e~ Xk Z 0 ) - 7 > z (qj) ^ 


\/ v Q q { z ) 

where D and are defined in Theorem [H 
Proof. We have 

(Qq,n (z k - e~ Xk Z 0 ) - 7Z 


F n 


(18) 


r 

n 

h 

V Q q J Z \ 


r 

n 

h 

V IJ Z \ 


n —1 


T {QqA Z ) ~7z(q)) + 




TRqJX.Z,) 


=0 


V Qq jZ) 


71—1 


+ 


E 


v l,S z ) 


v / ^KwI‘ =0 


RQ,q{\ Z k ). 


(19) 


oo 
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Hence 




dw 


dw 



( 

V Qq jZ) 



\ 



+ 

/ 

j n— 1 

- R Q,<1^ Z k) 

/ 

\ nE 

[ V U Z >] k -° 


L!(Q) 


Combining this with Q and (1141) we obtain (1161) . Similar argument leads to (1171) . 
Moreover, from (fT9l) . dTJ) and (ThO) we deduce (fISl) . ■ 


3.3 Power variation 

Consider, for every g £ N* even, the following power variation 

n— 1 

P q , n (Z):=-J2(Zi) q . ( 20 ) 

n i= o 

Define 

fe(«) := E[(Zo)’] = [^(Zo)] ,/2 - (21) 

We have the following almost sure convergence. 

Theorem 6 Suppose that Z is ergodic. Then, as n —>• oo 

P q ,n(Z) —► <$z(«) ( 22 ) 


almost surely. 

x 2 y 2 

Let 1): be the process defined in dT[) . Let c^fc = -p^yr j^oo e ^- x q H 2 k{x) dx be the coefficients 
of the monomial x q expanded in the basis of Hermite polynomials: 

q/2 

X q = ^C qj 2 k H 2 k {x). 
k =0 


10 






















Then we can write. 


V Pq jZ) = Vn{P q , n {Z) -E[(Z 0 ) q ]) 


[rzm q/2 _J 2 [ E 


[: rz(0 )] q/2 


1=0 \ 
n— 1 q/2 


z, 


(\Az(o) 


- E 


|_V\Az(o) / 


n 


EE% 2 fc(^)^ 2 fc(^i) 


i=0 fc=l 


fc=l V v i=0 


(g>2fc 


where Y) = Efe) = —M —. 

V 0) 

Furthermore, 


hL(z) 


9/2 (ok\\ n_1 

h(«)] , E<*( z )T : h Wi_i )l 

fc=l i,j=0 


2k 


9/2 


n—1 


2/c 


N( 0 )]’E 42 l=(Z)( 2 t)! 1 + -E( n - 1 -i)|rv(i)l 

/,-! \ ;-i 

9/2 / n—1 „ n—1 

b(0)l'E'h( Z )« ! ( l+2^|r,(jf --EjI’vWir ] ■ (23) 


k=l 


j =1 


1=1 


We will also need the following technical lemma. 


Lemma 7 Let (Zk)k ^0 a stationary Gaussian sequence with E(Zq) < 00 , and let A > 0 and 

q € N* even. Consider the sequence 

Rp, q { A, Z fc ) := (Z fc - e- Afc Zo)" - (Z fc ) 9 . 

Then for every p ^ 1 t/iere exits a constant c(A, ( 7 ) depending on A, (/ and E(Zq) such that 


Moreover for every e > 0 


almost surely as n —>• 00 . 


- n— 1 

- X] Rp ^' Z k) 

n ' 

fc =0 


ZX(Q) 


c(A,g) 

n 


1 

n £ 


n—1 

y, Rp,q(\ Zfe) 


/c=0 


0 


(24) 


(25) 
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Proof. The proof is straightforward by using similar arguments as in the proof of Lemma [4] ■ 
Applying Theorem [2] and Lemma [7] we conclude. 

Theorem 8 Let (Vp gn (^)) n>0 cind (f?p >9 (A, Z k )) n>0 the sequences defined in the above. 

1) Then there exist C depending on q and rz{ 0) such that 



On the other hand if l r r(^’)| 2 < °°> we can write 





(27) 


C 


vp q (Z) 


IvpAZ) 


n 


+ \ E 


V L 


(Y) jKfiVp 2 (Y)) + K 4 (Vp 2 (Y))+ v q (Z) — E 


vl 


XY) 


where vp q (Z) := lim ^,*,E 
2) Assume that ry(fc) = \k\ 


Vp qn (Z) . Moreover K±(V 2 ,n(Y)) verifies (EJ) and (EJ). 
“ with 0 < a < 7}. Then 


where D and are defined in Theorem [H 
Proof. Similar proof as in Theorem [8l ■ 


(28) 


3.4 Quadradic case 

In this subsection we suppose that q = 2. Then, in this case, we have 

^ n— 1 

Yq 2 JZ) = Vp 2 JZ) = J2{Zk-E[Zl]) 

Y n k =o 

= ^j^'Y J H 2 {Y k ) = E[Zl]V 2 ,fiY) 
= V 2 , n (Z). 

Thus, we obtain the following theorem. 

Theorem 9 We have 
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if \ r Y (^)| 2 < °o, there exists C > 0, ?iq ^ 1 such that for every n ^ no 


/ 


dry 


V 2 , n (Z) 


\ 


■;N 


E 


v,U z ) 


< c- 


(s 


|fc|<n K 


(E|fc|<nM*OF 


3/2 


n 


and hence also 


dw 


r 

n 

h 

v£ n { z) 


(Q2,n (Z k - e- Afc Z 0 ) - 7 z(2)) , N ) < 


C_ 

n 


( 29 ) 


2n 


!+ ( M*OI 3/2 

|fe|<n 


if\r Y (k)\ = \k\ 2 ; 


dry 


V 2 , n (Z) 


, N 


E 


ViUz) 


€ 


C 


(logn) 1 / 4 ’ 


and also 


dw 


n 


E 


vijz) 


^ Q 2 ,n (^Zfc 


-e~ Xk Zo) -7z(2)1 ,N I ^ 


C 


(logn) 1 / 4 ’ 


where in this case E [V^ 2 n (Z)] 4 h\ 2 {Z) log(n) = 4r|(0) log(n). 

if ry(k) = |A;| _ “ with 0 < a < -y, there exists C > 0 depending on a 


dry 


( 



,_ ^ 


v 2 ,n(Z) 

1 

/ 1 — 2a 

C 

< 

W £ 

V 2 2 n (Z) 

’ 2 \ 

D °° 

y 

■y/logn 


and also 


dw 


n 


E 


V 2 2 n (Z) 


(Q 2 ,n (z k -e- Xk Zo)- lz (2)),y 


(31) 


(32) 


(33) 


1 /l-2« \ C 

’2V -( 


where D and are defined in Theorem [H 
Moreover in this case E [V^ 2 n (Z)] ~ 4 ? (^ 2 2 ^ |r^(0)| 2 - 


is a 


Proof. Since — ! — — -' = " ! — =. then 

also 0, EU). Hence also, from (1191) . (12911 and (1 14(1 we deduce 


direct consequence of m Theorem 3] (see 


; (so) 
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Combining Q and © we obtain (13TT) . and also the estimates (1191) . (1311) and (1141) lead to (1321) . 
Now, suppose that ry(k) = \k\~ a with 0 < a < It is easy to see that (I23|) leads to 


E[V 2 2 n (Z)\ = 


r\j 


( n— 1 

l + 2j2\r Y (j)\ 2 


n 


1—2a 


1 - 2a 


rz( 0)| 2 . 




Thus, from |191 Theorem 5] and [11, Theorem 1], the estimate (I33|) is obtained. Hence also, from 
(fl9l) . (1331) and (fT4l) we deduce (l34l) ■ 


3.5 Improve the rate convergence 

Consider a centered stationary Gaussian process Z = (Zk) kGZ with covariance rz{k ) = E^Z^Zk) < 
oo and \rz(k)\ ^ 1 for k € Z. 

Define the centered stationary Gaussian process Z ^ as follows: 

z\, ^ = Zk+i — Z^, k £ Z 

and for every p ^ 2 

4 P) = zfc 1 ’ - zC'D € Z. 

Now, applying Theorem [2] we conclude the following result. 


Theorem 10 Assume that k a rz(k ) converges to a constant for some a < — Then for every 
p ^ 2 there exists C depending on p , q and rz{ 0) such that 


ITV 


V Qq , n ((Z W)) 


-;N 


E 




< Cn a ~ 2p+ 2 . 


(35) 




In particular, if a = 2 H — 2 f/ien /or every p^ 2 there exists C depending on p, q and rz( 0) such 
that 


( 


drv 


V Qq J(Z (p) )) 


\ 


,N 


E 


V IJ( Z{P) )) 


Cn 2H-2 P -i' 


(36) 


This leads that if p ^ 2, then 




is asymptotically normal for every H € (0,1). 
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4 


Applications to Ornstein-Uhlenbeck processes 

4.1 Fractional Ornstein-Uhlenbeck process 

In this section suppose that X = {Xf,t ^ 0} is an Ornstein-Uhlenbeck process driven by a frac¬ 
tional Brownian motion B H = {77^,7 ^ 0} of Hurst index H G (0,1). That is, X is the solution 
of the following linear stochastic differential equation 

X 0 = 0; dX t = —6X t dt + dBf 1 , t ^ 0, (37) 

where whereas 6 > 0 is considered as unknown parameter. 

The solution X has the following explicit expression: 

X t = [ e-W-^dB?. 

Jo 

We can also write 


X t = Z? - e~ dt Z? 


(38) 


where 


z e t = 


f 

J —c 


e~ e{t - s) dB 


H 


Moreover, Y e is an ergodic stationary Gaussian process. 

We will need the following result. 

Lemma 11 Let H G (0, ^) U (^, 1], m,m' > 0 and —oo ^a<b^c<d< oo. Then 

E ( ( b e ms dB H (s) [ d e m,t dB H {t) \ =77(217-1) f dse ms f dte m '\t - s) 2H ~ 2 
\J a Jc J J a Jc 

Proof. We use the same argument as in the proof of |T0l Lemma 2.1] ■ 


Lemma 12 Let H G (0, i) U (^, 1), m,m! > 0 and let Z e be the process defined in \39\) . Then, 


and for large |t| 



= HT(2H)d~ 2H 


E 


r/Q ry6 


77(277 - 1) 

e 2 1 


\2H-2 


Proof. see ( TO] Theorem 2.3] or Lemma [171 ■ 


(39) 


(40) 
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4.2 Construction and asymptotic behavior of the estimators 

Fix q ^ 2 and assume that q is even. 

From (I38[) we can write 

1 n— 1 

Q q ,n( X ) = Q q ,n(Z 9 ) + ~J2 R Q^^ Z k) ( 41 ) 

1 k =0 

Combining (1411) . Lemma [4] and the fact that Z e is ergodic we conclude that, almost surely, 

lirn Q q ,n{X) = lirn Q qn (Z e ) 

n — >oo n —>-oo 

= 7 zo(q) 

= (psT? - 1 ) ,/2 

:= Hq(9). 

Hence we obtain the following estimator for 9 

0 q ,n = 9-q 1 [Q q ,n(X)] . 

As consequence, we have the following strong consistence of 9 q ^ n . 


Theorem 13 Let H G (0,1). Then, as n —> oo 

0 q , n —► 9 (42) 

almost surely. 

Combining (1411) and Theorem [8] we conclude the following result. 


Theorem 14 Denote N ~ A/"(0,1). If H G (0, |], then there exists C depending on q, H and 9 
such that 


d w 


n 


E Kn{Z d )\ 


Tq Kn =7 Cl 


In particular, 


if He (o, f) 


_i 5 

n 4 , if 0 < H < - 
8 

n - 4log4(n), if H = — 

8 

9 TT _ ^ 5 3 

ra 2, */ - < H < - 

log ~H n ), i f H = \- 


n[lig(9 q , n ) -f* q (0)) ^ AC (0,<r^(Z f 


2/,70 


(43) 


( 44 ) 


where a 2 (Z e ) = b 2 g>2k (Z d )(2k)l (l + 2 = lim, woo E [V 2 n (Z e )], 
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ifH = | 


(/% (e„n) - (.,(»)) ^ jv (o,< 2 (z»)) 


where in this case E [V^ n (Z 0 )] ~ 462 2 (^)log(n). 
7n the case when H € (|, 1), we have 




1 


n 2 

where is defined in Theorem [H 

Thus we deduce the asymptotic distribution of 9 qyn . 

Theorem 15 If H € (0, |], then 


Vd 


n 


E [V£ n (Z°j\ 


e q , n -e)^M(o,(n' q (6)) 2 ) 


If H € (|, 1), fherr 


Proof. We can write 


n 2H ~^ 


3 l 


_ ^ j kg,2 ^ 






n q,n) hq — /hj(£g,n)V^ {®q,n 9^j 


where is a random variable between 0 and 

Combining this with Theorem 1141 we obtain the desired conclusion. ■ 

Quadratic case: Now we will discuss the particular case when q = 2. 
Define 


If := 






z? 


\FwV) ^[(Z 0 e )2] fW 

Combining (141 [) and Theorem [9] we conclude the following result. 
Theorem 16 If H € (0, |] ; £hen 

1 


,-2H 


/ 




HT(2H) 


\ 


E 


vU z9 ) 


■V^(of 2 n H -e- 2H ),N 


< C x < 


VS- 

‘ff !/ "=3 


b ,H T 

_ l _ if H = -. 

, (log(n)) 1 / 4 ’ 4 


(45) 


(46) 


(47) 


(48) 


(49) 


(50) 
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As consequence, for every H € (0, |) 


n 02,n - O 


AMO 


and if H = 


n 


log(n) 


'2,n 


0 iH+2 a 2 (Z e ) 
(HT(2H + l)) 2 




where in this case E [VjP n (Z 0 )] ~ (r(277 + 1)0 2 ^) log(n) 

// § < if < 1, 


(51) 


(52) 




i?r(2i?) 






where E [V£ n (Z e )\ ~ w 
As consequence 


Vln(Z G ) 

(r(2i7 + l) 0 " 2i/ ) 2 . 


’ 2 


D 




C 


\/log(n) 


V^(^-02,n) ^ ^oo. 


4.3 OU driven by fractional Ornstein-Uhlenbeck process 


(53) 


(54) 


In this section suppose that V = {A^, t ^ 0} is an Ornstein-Uhlenbeck process driven by fractional 
Ornstein-Uhlenbeck process V = {14, f ^ 0} given by the following linear stochastic differential 
equations 

( X 0 = 0; dX t = —9X t dt + dV t , t > 0 
\ Vo = 0; dV t = -pV t dt + dBf 1 , t ^ 0, 

where B H = {B^,t ^ 0} is a fractional Brownian motion of Hurst index H € (0,1), whereas 
0 > 0 and p > 0 are considered as unknown parameters such that 6 ^ p. 

For m we have the following results related to Xp. 


X, = 


p xf + J—x? 


p-q 0-p * 


where for m > 0 

X™= f e- m(t - s \lB^. 

Jo 

On the other hand, we can also write the system (1551) as follows 


(56) 


(57) 


dX t = — (9 + p) X t dt - p9E t dt + dB 


H 
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where for 0 ^ t ^ T 


l 


4 v t - x t _ X f - x p t 


= / X,ds = 


P-0 


(58) 


We can also write 


-ym ryrn „—mt rym 

X t - - e 


(59) 


where 


Z\ n = j 


(60) 


As consequence 


W = 


p „ n 6 „ fl / pe pt „„ #e et 


Z p + ^—Z e f - ( Z, e 


p-d * d-p 1 \p — 0 u 9 — p 

6 


■ = Z e t ’ p +-^-R{p,ZP) + -^—R{6,Z°) 

p — u u — p 


(61) 


and 


S* = 


On the other hand, the process 


Z e t - Z p t e~ 9t Z° - e-<*Z? 

p — 9 p — 6 


:= E 




(62) 


^rri 


(63) 


is an ergodic stationary Gaussian process. 


Lemma 17 Let H E (0, 2 ) U (g, 1), vn,m' > 0 and /et Z m be the process defined in \60\) . 


A (m, m!) := E 


rym rym 
^0 


HT(2H) ( x _ 


m + m‘ 


) f 1-2 H I ( /\1-2H\ 

j ym + (m ) J 


Then, 

(64) 


and /or /arge |t| 


E 


rym rym 
^0 


H(2H — 1) U| 2if_2 


mm 




(65) 


This implies that for H € (0, g) U (^, 1) 


Vx(0,p) := E 


7&, P 


HT(2H) 2 _ 2H n2-2H] 


p 2 — 6 2 


■\ P 2 ~ 2H - e 2 - 2M ], 


rm(0,p) -=E 


^0 


HT(2H) 2 H -2 H-I 

- p 
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and 


E 


y6,p-^6,p 

^0 


= 0. 


Proof. By using |10l Proposition A.l], we can write 

r 0 p 0 


rym rym 
^0 ^0 


= mm 

mm‘ 


2 

mm‘ 


' f f e mu e rn ' v E (B^Bjf) dudv 
J —oo J —oo 

e mu e m ' v {{-u) 2H + {-v) 2H -\ V - U \ 2H ) dudv 


— OO J —OO 

/ r 0 r 0 


—oo ./ —oo 
/ /»oo /*oo 


/*oo /*c 

Jo Jo 


e e 


{u 2H + 


v — \v — u\ 


r ) dudv 


F(2H + 1) f 1 _ 2H ( ,^i-2 H\ 

— -— \ m +77i 

2(m + m') V / 


Thus the estimate (1M1) is proved. 
Now, let 0 < e < 1 


E ( Z^Z™ 


= e 


(/" 

r 

\J —00 

J—1 

(/" 

e mu dB^ d 

\ J 00 

J — 1 


e m ' v dB 7 


= e 
:= A + B 

where, using ] T0l Proposition A.l] it is easy to see that 

\A\=0 (e"” 1 '*) . 

On the other hand, by Lemma [IT] and integrations by part 


rt 
f et 


e mu dB% / e m ' v dBf 


B = H{2H - l)e~ m,t j du e mu f dv e m ' v {v - u) 2H ~ 2 

J —oo Jet 

rO rt—u 

= H(2H - l)e~ m,t / du e mu / dz e TO '( u +*) z 2H ~ 2 


f —OO 
roo 


' et—u 


roo rOA(t—z) 

H(2H - l)e- mlt / dz e m ' z z 2H ~ 2 / du 7 m+m > 

Jet J et—z 

H(2H-V)_{ t°° e -m{z-t) z 2H-2 dz+ f* e -m'(t-z) z 2H-2 dz + e -m't{l-s) f°° e ~m(z-et) z 2H-2 dz 

\Jt Jet Jet 


m + m 

H(2H - 1) ft 2H ~ 2 2H-2 


(m + in') \ 777 
(et) 2H ~ 2 


+ 


e -m(z-t) z 2H-3 dz + 


2H-2 


m 


m 


\-\2H-2 9 7 _r _ 9 rt r oo \ 

_ e -m'(l-e)t _/ e -m'(t-z) z 2H-S ( p z e -m't(l-e) / e ~m(z-et) z 2H-2 j 

m’ rri J £t J £t ~ 7 


g(2g ~ l)7 g - 2 + o (t 2 ^- 2 ) , 

777.777/ v 7 


20 

















the last inequality comes from the fact that 


poo poo 

/ z 2H ~ 3 dz < r 1 / e~ my dy 

Jt Jo 

—y 0, as t —>• oo, 


and 


t 2-2H ^ e - m '(t-z) z 2H-3 dz ^ £ 2H-3 r l / e ~m'(t-z) dz 

J et J et 

-(l-e)i 


ft 


'7 


= e 2 ^- 3 *" 1 / e- m ' y dy 


0, as t —> 00 , 


t 2-2H e -m’t(l-e) ^ 0? as t ^ OO. 

So, we conclude that the estimate (165[) is obtained. ■ 


4.4 Construction and asymptotic behavior of the estimators 
Case: X and S are observed 

Combining (l6lj) . (l62j) . (ITTI) . (l63|) and Lemma 0] we conclude that 

lim (Q q , n (X),Q q>n (£)) = lim (Q 9 , n X 0 ’ P ), QqA^ P )) 

= (TV®.#- (9), («)) 

= , g J 2q/2 (7v(0, p) - l ) q/2 , (»7e( 0, p) - 1) 9/2 ) 

:= 4(M- 

Hence, we obtain the following estimator for (0,p). 

{^q,m Pq,n ^ — G q (Qq,n(-^) i Qq,n(^)) . 
where G q is the inverse function of F q . 

By construction, we have the following strong consistence of | 'o q>n , p q . 

Theorem 18 Let H G (0,1). Then, as n —> 00 

(0q,m Pq,n^ t iP i P') 

almost surely. 

Combining m, (ed, Lemma DU and Theorem [8] we conclude the following result 


( 66 ) 
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( 67 ) 


Theorem 19 Let H € (0, |). Then 

[(Q,,nW, Q,,n(S)) - (0, p)\ A A7 (0, r(0, p)) 


where 


where 


mp) = 


a q(Y e,p ) a q (Y^,E°’P) 
a q {Y e ’P,Y0,p) a 2 ( E^) 


( 68 ) 


a 2 q (Y e ’P) = Var [H q (Y°’ p j\ + 2 £ Cov [H q (Y°’ p ), H q (Y°' p ) 

k =1 
00 

a q(^ e,P ) = Var [H q (^ p )] +2^ Cov [h^), H q (Y^ p ) 

k= 1 

^ (Y^,sM = Cov \H q {Y°'P),H q {Y e *) 1 + ]T Cot; [ff,(y o ^),^(E^) 


fcez* 


Theorem 20 Lei iL € (0, |). TTten 


9q,n- d ’P q ,n- p) ■A/' (0, Jc q {vx, Vs) T(9,p) 1 Jc q {vx{9, p), Vs{9, p))) (69) 


Proof. By Taylor’s formula we can write 


H \^@q,n @Pq,n Pj — \Z^ i_Qq,n{_-^-) Tf^’P ((?)> Qq,n(^) Te {.VX i ??e) H" d n 

where d n converges in distribution to zero, because 

IKII < c(0,p,L7)v/i;||Q ?>n (X) - 7 ys, P (g),Q 9in (S) - 7s 0 .e(9))H 2 - 
It is easy to see that if for any w S f l there exists uq(iv) € N such that X n (w) = Y n (w) for all 

n ^ no(w ) and 0 as n —»• oo, then T n —> 0 as n —»• oo. 

Combining this and above estimates we obtain the desired result. ■ 

Quadratic case: Here we suppose that q = 2. In this case we have: 

L 2 is a positive function of the variables (x, y ) in (0, + 00) 2 defined by: for every (x, y ) € (0, + 00) 2 

vrry? (y 2 ~ 2H ~ x 2 ~ 2H , x~ 2H — y~ 2H ) if x 7 ^ y 


F 2 (x,y) =HT(2H) x f 


((1 — H)x 2H ,Hx 2H 2 ) if y = x. 


(70) 


Since for every (x, y) € (0, + 00) 2 with x / y 
Jf 2 (x,y) = T(2 H + 1) 


' (l—H)x 1 ~ 2H (x 2 —y 2 ')—x(x 2 ~ 2H —y 2 - 211 ') {l-H)y'~ 2H (y 2 -x 2 )-y(y 2 ~ 2H -x 2 ~ 2H )' 

- 2 / 2^2 ^^>2 _ y 2 , y 2 

i/ai- 2ff - 1 (ai 2 -y 2 )+ir(ai- 2H -i / - 2 ^) Hy-™- 1 _ X -2H ) 

{x 2 -y 2 ) 2 (x 2 -y 2 ) 2 

the determinant of Jf 2 (x,y) is non-zero on in (0,+oo) 2 . So, F 2 is a diffeomorphism in (0,+ 00) 2 
and its inverse G 2 has a Jacobian 

, Hy- 2 »-i(y 2 -x 2 )+y{y- 2H -x- 2 ") (1 - H )y'~ 2H (y 2 -x 2 )-y{y 2 - 2H -x 2 ~ 2H )' 

r( 2 iL + 1 ) I (x 2 -y 2 ) 2 (x 2 -y 2 ) 2 


Jg 2 («> b ) = 


det Jp 2 (X, y) 1 _ Hx- 2 »-yx 2 -y 2 )+x(x- 2H -y- 2 ») (1 -H)x^~ 2H {x 2 -y 2 )-x(x 2 ~ 2H -y 2 ~ 2H ) 

{x 2 -y 2 f (x 2 -y 2 ) 2 


where (x, y) = G 2 (a, b). 
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4.5 Fractional Ornstein-Uhlenbeck process with the second kind 

Let U = {Ut, t ^ 0} be a fractional Ornstein-Uhlenbeck process with the second kind defined as 

Uo = 0, and dUt = —aUtdt + dY^ l \ t ^ 0, (71) 

where Y^ = Jq e~ s dB as with a s = He h and B = {Bt,t ^ 0} is a fractional Brownian motion 
with Hurst parameter H £ (^, 1), and where a > 0 is a unknown real parameter. 

The equation (1711) admits an explicit solution 

U t = e~ at [ e as dYp = e~ at [ e^~^ s dB as 

Jo Jo 

rat 

= H^-^ H e~ at / r {a ~ 1)H dB r . 

J ao 

Hence we can also write 


U t = U? + R(a, U a ) 


(72) 


where 


U? 



e^ a - l >dB aa 


fj(l-a)H e -ott 



A a ~ 1)H dB r . 


Lemma 21 Let H £ (4 1). Then, 


E 


(US 


(2 H - 1 )H 2H 


a 


P(l-H + aH, 277-1). 


and for large |t| 

ru°(t) = E [t/ 0 Q U“] = O (e- rnin{Q ^ }/ ) 
Proof. We prove the first point (1731) . We have 


E 


(US 


) 2 ] = 77(277 - [ ° dyy (a - ])H [ ° dx x^~ 1)H \x - y\ 2H ~ 2 

- 1 Jo Jo 


rao 


= 277(277 — l)77 2 ( 1 ” Q )^ / dyy (a ~ 1)H dx x ^ a ~ l)H {y - x) 


2H-2 


rao H 

277(277 - \)H 2 U~ a ) H / dyy 2aH ~ l / dz (l - z) 

Jo Jo 

(277 - 1)H 2H 


2H-2 


a 


p(l- H + aH, 277- 1). 


(73) 


(74) 


Thus (1731) is obtained. 

For the point (174)) see [16]. 
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Combining (|72[) . Lemma [4] and the fact that U a is ergodic we conclude that almost surely 


lim Q q: n(U) = lim Q q , n (U a ) 


n —>-00 


n—>00 


= iu<* ( q ) 

q\ 


(§)! 2<?/2 


;|)!29/2 


E 


(u§ 


\?/2 

-0 




a 


■ = Vg(a). 

Hence we obtain the following estimator for a 

Otq,n = Vg \Qq,n(U )] . 

By construction, we have the following strong consistence of a qtT 
Theorem 22 Let H G ( 5 , 1 ). Then, 


as n —> 00 


a , 




a 


(75) 


almost surely. 

Combining fl72|) and Theorem [HI Lemma [21] and Lemma [4] we conclude the following result 

Theorem 23 Denote N ~ jV(0, 1). If H G ( 5 , 1 ), t/ien t/iere exist C depending on q, H and a 
such that 


dw 


n 


E [V 2 n (U«)\ 


{Vq (a q , n ) - v q (a )), N ^ Cn 4 . 


In particular, 


\fn(v q (a qin ) — v q (a)) —> (0, a 2 (U a )) 

TdPA _ m 


(76) 


(77) 


where a 2 (U a ) = £^1 # q ,2 k (U a )(2k)\ (l + 2 ££1 ) = lim™ E [V 2 n (U a )\. 


«=! 

Theorem 24 Let fL G (~, 1). T/ien 


/£ wiwi\ (S, - n ~ a) ^ M (°- M ( “)) 


-2 


(78) 


Proof. We can write 


y/n (v q (a q ,n) ~ *4/(a)) = ^Mq,n)y/n (S 9i „ - a) 


where £ 9jri is a random variable between a and S 9jn . 

Combining this with Theorem 1231 we obtain the desired conclusion. ■ 

Quadradic case: Here we suppose that q = 2. We will study the asymptotic distribution of 

®2,n- 
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Theorem 25 Let H G (i, l). Then there exists C depending on H and a such that 



n 


C 


(y 2 (a 2 ,n) - v 2 (ol)) , N ^ — 

vl „(£/«) I 


( 79 ) 


Proof. From Lemma [2T1 we have for any H € (^, 1), \ r Z a {k)\ 3 / 2 < oo and YlkeZ \ r Z a {k )\ 2 < 

oo. Hence by Theorem O the conclusion is obtained. ■ 
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